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f-*) \ The cranked relativistic Hartree+Bogoliubov theory has been apphed for a 

' systematic study of the nuclei around ^^^No, the heaviest elements for which 

■ detailed spectroscopic data are available. The deformation, rotational re- 

Cf^ , sponse, pairing correlations, quasi-particle and other properties of these nuclei 



have been studied with different parametrizations for the effective mean-field 
Lagrangian. Pairing correlations are taken into account by a finite range 
two-body force of Gogny type. While the deformation properties are well 
reproduced, the calculations reveal some deficiencies of the effective forces 
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' both in the particle-hole and particle-particle channels. For the first time, the 

quasi-particle spectra of odd deformed nuclei have been calculated in a fully 
\ self-consistent way within the framework of the relativistic mean field (RMF) 

;h ' theory. The energies of the spherical subshells, from which active deformed 

^ states of these nuclei emerge, are described with an accuracy better than 0.5 

MeV for most of the subshells with the NLl and NL3 parametrizations. How- 
ever, for a few subshells the discrepancies reach 0.7-1.0 MeV. In very heavy 
systems, where the level density is high, this level of accuracy is not suffi- 
cient for reliable predictions of the location of relatively small deformed shell 
gaps. The calculated moments of inertia reveal only small sensitivity to the 
RMF parametrization and, thus, to differences in the single-particle struc- 
ture. However, in contrast to lighter systems, it is necessary to decrease the 
strength of the DIS Gogny force in the pairing channel in order to reproduce 
the moments of inertia. 

PACS: 21.60.Cs, 21.60Jz, 27.90.+b, 21.10.Pc 
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I. INTRODUCTION 



The possible existence of shell-stabilized superheavy nuclei, predicted with realistic nu- 
clear potentials [1-3] and the macroscopic- microscopic (MM) method [4-6] , has been a driv- 
ing force behind experimental and theoretical efforts to investigate the superheavy nuclei. 
These investigations pose a number of experimental and theoretical challenges. The recent 
discovery of elements with Z=112 [7], Z=114 [8] and Z=116 [9] (for review of the present ex- 
perimental situation see Refs. [10-12]) clearly shows great progress on the experimental side, 
but also indicates difficulties in the investigation of nuclei with low production cross-sections 
and analyses based only on 1 or 2 events. 

The theoretical challenges are also considerable since different theoretical methods pre- 
dict different spherical shell closures. Modern calculations based on the MM method with 
the Woods-Saxon [4,13,14], Nilsson [5], and folded Yukawa [15] potentials indicate Z = 114 
and = 184 as the spherical shell closures. It is necessary to say, however, that some earlier 
calculations indicated Z — 126 as a possible magic number (see Ref. [10] for a review). There 
are differences in the predictions of self-consistent calculations, which depend both on the 
approach and on the effective force. Self-consistent calculations based on the Hartrce-Fock 
method with Skyrme forces (SHF) predict spherical shell closures at Z = 126 and N = 184 
for most of the forces [14,16,17]. However, some forces indicate Z — 114 (SkI4) and Z — 120 
(SkI3) as proton shell closures, while some predict no doubly magic superheavy nuclei at all. 
On the other hand, the rclativistic mean field theory (RMF) prefers Z — 120 and N = 172 
as spherical shell closures [16,17]. However, Z = llA and = 184 also appear as the shell 
closures in some RMF calculations [18,19]. For a detailed comparison of the predictions 
of the different Skyrme and RMF calculations, see Refs. [17,20]. Hartree-Fock-Bogohubov 
(HFB) calculations with the Gogny force give Z — 120, 126 and N — 172, 184 as spherical 
shell closures [21]. 

Clearly, the accuracy of predictions of spherical shell closures depends sensitively on the 
accuracy of describing the single-particle energies, which becomes especially important for 
superheavy nuclei, where the level density is very high. Variations in single-particle energy 
of 1 — 1.5 MeV yield spherical shell gaps at different particle numbers, which restricts the 
reliability in extrapolating to the unknown region. 

Usually, the MM method describes the single-particle energies rather well. This is due 
to the fact that the experimental data on single-particle states arc used directly in the 
parametrization of the single-particle potential. Moreover, different paramctrizations of 
the single-particle potential are used in different mass regions. However, the extrapolation 
of the single-particle potential may be much less reliable since it is not determined self- 
consistently. For example, microscopic models predict that the appearance of the shell 
closures in superheavy nuclei is influenced by a central depression of the nuclear density 
distribution [17,22]. This effect cannot be treated in self-consistent way in current MM 
models. 

Although the nucleonic potential is defined in SHF and RMF approaches in a fully self- 
consistent way, this does not guarantee that single-particle degrees of freedom are accurately 
described, as indicated by the large variety of the paramctrizations (more than 60 for SHF 
and about 20 for RMF [20]). In addition, the parameters have been fitted in almost all cases 
to the bulk properties of spherical nuclei. Single-particle information on spin-orbit splittings 
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is used only in the fits of the parameters of the Skyrme and Gogny forces. The spin-orbit 
interaction is a relativistic effect, which arises naturally in the RMF theory. Thus, available 
RMF fits were obtained without the use of any single-particle information. For heavy nuclei, 
the calculated RMF single-particle states were directly compared with experimental data 
only in spherical nuclei (see, for example, Refs. [23,24] and references quoted therein). These 
comparisons, however, do not reveal the accuracy of the description of the single-particle 
states because the particle- vibration coupling, which can affect considerably the energies of 
single-particle states in spherical nuclei [25-27], has been neglected. 

Compared with the MM method, self-consistent calculations have been confronted with 
experiment to a lesser degree and for a smaller number of physical observables (mainly 
binding energies and quantities related to their derivatives). For many parametrizations, 
even the reliability of describing conventional nuclei is poorly known. In such a situation, 
it is important to perform a comprehensive study of the heaviest nuclei for which detailed 
spectroscopic information is available. The results of such a study will allow us to better judge 
the reliability of predictions for superheavy nuclei. The experimental data on deformed nuclei 
around ^^^No provides sufficient information for such a test. The purpose of this work is to 
compare the predictions of RMF theory with these data. 

RMF calculations have been compared with experiment in this mass region in only Ref. 
[28]. However, the comparison was restricted to binding energies and quantities related to 
their derivatives. In addition, the pair correlations were treated in the BCS approximation, 
with no particle number projection. 

In the present manuscript, the cranked relativistic Hartree-Bogoliubov (CRHB) theory 
[29,30], with approximate particle number projection by means of the Lipkin-Nogami method 
(CRHB-I-LN theory), is employed for a detailed investigation of a wide set of experimental 
observables. The use of the Lipkin-Nogami method has the clear advantage of avoiding the 
collapse of pairing correlations at large shell gaps. We address for the first time the question 
of the blocking procedure in odd mass nuclei in the framework of the RMF theory, with 
effects of time-reversal symmetry breaking taken into account in a fully self-consistent way. 
The calculated binding energies, deformations, moments of inertia and quasi-particle states 
are compared with experiment. 

The paper is organized in the following way. In Sect. II, a brief description of the 
CRHB-I-LN theory and of some specific features of the present calculations are given. In 
order to outline the general features of the evolution of physical observables as a function of 
proton and neutron number, systematic calculations with different RMF parametrizations 
are performed along the Z — 100 (Fm) isotope and N — 152 isotone chains. In addition, 
calculations are carried out for Cm, Cf and No nuclei, for which experimental data are 
available. The rotational properties are studied in detail in Sect. HI, the deformations are 
discussed in Sect. IV and the shell structure in Sect. V. The quasiparticle spectra of selected 
nuclei are compared with experiment in Sect. VI. Finally, Sect. VII summarizes our main 
conclusions. 

II. THEORETICAL FORMALISM 

In relativistic mean field (RMF) theory [23,31,32] the nucleus is described as a system of 
point-like nucleons (Dirac spinors) coupled to mesons and photons. The nucleons interact 
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by exchanging scalar a-mcsons, vector uj-, p- mesons and photons. The isoscalar-scalar a- 
mesons generate strong intermediate range attraction between the nucleons. For the vector 
particles we have to distinguish the time-like and the space-like components. In the case 
of photons, they correspond to the Coulomb field and the magnetic field when currents 
are present. For the isoscalar-vector a;-meson, the time-like component provides a very 
strong repulsion at short distances for all combinations of nucleons, pp, nn and pn. For 
the isovector-vector p-meson, the time-like components give rise to a short range repulsion 
for like nucleons {pp and nn) and a short range attraction for unlike nucleons (np). They 
also have a strong infiuence on the symmetry energy. The space-hke components of the u- 
and p-mesons lead to an interaction between currents, which is attractive in the case of the 
a;-meson for all combinations {pp, nn and pn) and in the case of the p-meson attractive for 
pp and nn, but repulsive for pn. Within the mean field theory, these currents only occur in 
cases of time-reversal breaking mean fields, which applies to rotating or odd-mass nuclei. 

The cranked relativistic Hartree-Bogohubov (CRHB) theory [29,30] extends the RMF 
theory to rotating nuclei and includes pairing correlations. If an approximate particle num- 
ber projection is performed by means of the Lipkin-Nogami (LN) method [33-36], the ab- 
breviation CRHB+LN will be used. Since the theory is described in detail in Ref. [30], only 
the features important for the present discussion will be outlined below. 



A. The CRHB-I-LN equations 

The CRHB-I-LN equations for the fermions in the rotating frame are given by [30] 

(h'^-X'-nJ, A \fU{r)\ _j,,fU{r)\ 

[ -A* -h'^ + X' + nj;)[v{r)), '''[V{r)),^ 

where 

h'i,^hD + ^X2p-2X2Tr{p), (2) 
A' = Ai + 2A2 , (3) 
Ei^Ek-X^. (4) 

Here, Hd is the Dirac Hamiltonian for the nucleon with mass m; Xi is defined from the 
average particle number constraints for protons and neutrons; Pr = V*V^ is the density ma- 
trix; Uk{r) and T4(r) are quasiparticle Dirac spinors; E]. denotes the quasiparticle energies; 
and Jx is the angular momentum component. The LN method corresponds to a restricted 
variation of A2((AA^)^) (see Ref. [30] for definitions of Ai and A2), where A2 is calculated 
self-consistently in each step of the iteration. The form of the CRHB-I-LN equations given 
above corresponds to the shift of the LN modification into the particle-hole channel. 
The Dirac Hamiltonian ho contains an attractive scalar potential S{r) 

S{r)=g,a{r), (5) 

a repulsive vector potential V^)(r) 

1 - rs 

Vo{r) = g^ujQ{r) + gpT3po{r) + e ^ Ao{r), (6) 
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and a magnetic potential V{r) 



V{r) = g^oj{r)+gpnp{r) + e^-^A{r). (7) 

The last term breaks time-reversal symmetry and induces currents. In rotating nuclei, the 
time-reversal symmetry is broken by the Coriolis field. Without rotation, it is broken when 
the time-reversal orbitals are not occupied pairwise. In the Dirac eqiiation, the space-like 
components of the vector mesons uj{r) and p{r) have the same structure as the space- 
like component A[r) generated by the photons. Since A(r) is the vector potential of the 
magnetic field, by analogy the effect due to presence of the vector field V{r) is called 
nuclear magnetism [37]. It has considerable infiuence on the magnetic moments [38] and 
the moments of inertia [39-41]. In the present calculations the spatial components of the 
vector mesons are properly taken into account in a fully self-consistent way. The detailed 
description of the mesonic degrees of freedom in the CRHB-I-LN theory is presented in Ref. 
[30]. 

The pair field A is given by 

A^A^, ^ ^-Y^VZa^c (8) 

cd 

where the indices a,b, . . . denote quantum numbers which specify the single-particle states 
with the space coordinates r, as well as the Dirac and isospin indices s and r. It contains 
the pairing tensor k 

K^V*U^ (9) 

and the matrix elements of the effective interaction in the particle-particle {pp) channel, 
for which the phenomenological non-relativistic Gogny-type finite range interaction 

VPP{1, 2) = / ^ e-KTi-ra)//..]' x {Wi + BiP" - HiP^ - M^P^P^) (10) 

i=l,2 

is used. The clear advantage of such a force is that it provides an automatic cutoff of high- 
momentum components. The motivation for such an approach to the description of pairing 
is given in Ref. [30,42]. In Eq. (10), W,, Sj, and Mj (i = 1,2) are the parameters 
of the force and P'^ and P'^ are the exchange operators for the spin and isospin variables, 
respectively. Note that a scaling factor / is introduced in Eq. (10). In our previous studies, 
the original (scaling factor / = 1.0) parameter set DIS [43,44] provided a good description 
of the moments of inertia in the A ~ 75 [45], A ~ 160 - 170 [46] and A ~ 190 [29,30] mass 
regions. As discussed in Sect. Ill B, it produces pairing correlations in the A ~ 250 mass 
region that are too strong, and, thus, it has to be attenuated (/ < 1.0). 

As a measure for the size of the pairing correlations in Hartree-(Fock)-Bogoliubov calcu- 
lations, we use the pairing energy 

Epairing = " ^Tr(AK). (11) 
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B. The RMF parametrizations 



In the present study, the NLl [47], NL-Z [48], NL3 [49], NLSH [50] and NL-RAl [19] 
parametrizations will be compared in order to see how well observables, such as the mo- 
ments of inertia, the deformations, the quasiparticle energies, the separation energy and the 
quantity S2n{Z,N) related to its derivative, agree with each other and with experiment. 

These sets differ in the experimental input used in the fitting procedure. The binding 
energies of a number of spherical nuclei were included in the fit of all those sets, but the 
selection of nuclei was different. NLl and NL-Z employ the data mainly from the valley of 
beta-stability, while additional information on neutron-rich nuclei has been used in the fit of 
the NL3 set. Moreover, there is a difference in the selection of additional observables used 
in the fit. Charge diffraction radii and surface thicknesses were included in the fit of NLl 
and NL-Z sets [47,48]. The NL-Z set is a re- fit of NLl where the correction for spurious 
center-of-mass motion is calculated from an actual many-body wave function [48]. On the 
contrary, NL3 and NLSH employ data on charge and neutron radii [49,50]. This (together 
with the fact that in the NL3 set more experimental data on neutron rich nuclei were used in 
fitting procedure) provides better a description of isospin, surface and symmetry properties 
of finite nuclei in the NL3 and NLSH sets. Unfortunately, Ref. [19] does not state which 
data the NL-RAl set is fitted to. 

The sets NLl, NL3 and NLSH have been used extensively in RMF studies and tested 
on a wide range of physical observables related, for example, to the ground state properties, 
rotational properties, properties of giant resonances etc.; see Ref. [23] for review. The sets 
NL-Z and NL-RAl have been tested only for observables related to ground state properties. 
The set NL-Z is a re-fit of NLl with a correction for spurious center-of-mass motion [48] 
given by 

Ecm = -{Plj/2mA, (12) 

where Pcm is the total momentum operator in the center-of-mass frame, m is the nucleon 
mass and A the mass number. This term is added after the variation is performed to cir- 
cumvent two-body terms in the mean field equations. Thus, the use of other prescriptions 
instead of Eq. (12) for the treatment of the center-of-mass motion with NL-Z will affect 
only the binding energies and the quantities related to their derivatives. In all our calcula- 
tions (including those with NL-Z), the correction for the spurious center-of-mass motion is 
approximated by its value in a nonrelativistic harmonic oscillator potential 

E,^^~A1A-'/^ MeV. (13) 

This is consistent with the NLl, NL3, NLSH and NL-RAl parametrizations. As illustrated 
in Ref. [51], Eq. (13) is a very good approximation to Eq. (12) in the A ~ 250 mass region: 
the difference between two prescriptions does not exceed 0.3 MeV, which is only RiO.017% 
correction to the typical binding energy and changes smoothly with the mass number A. 
Based on the results given in Fig. 2 of Ref. [51] one can estimate that in this mass region 
the use of Eq. (13) instead of Eq. (12) will affect two-particle separation energies S2n{Z, N) 
and S2n{Z, N) by at most 0.030 MeV. This justifies the use of Eq. (13) for NL-Z. 
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The parametrization NL-RAl has been introduced recently in Ref. [19]. A number of 
conclusions of this article has been strongly questioned in Ref. [52] . in part due to the use of 
an unrealistically strong pairing interaction. However, if a more realistic pairing is employed, 
this parametrization provides a rather good description of the binding energies (see Fig. 2 
in Ref. [52]). 

C. Details of the calculations 

The CRHB(+LN) equations are solved in the basis of an anisotropic three-dimensional 
harmonic oscillator in Cartesian coordinates with the deformation parameters /?o = 0.3, 
7 = 0° and oscillator frequency hcuo = 41A~^/^ MeV. All fermionic and bosonic states 
belonging to the shells up to A^^ — 14 and Nb — 16 are taken into account in the diag- 
onalization of the Dirac equation and the matrix inversion of the Klein-Gordon equations, 
respectively. The detailed investigation of 246,248,250pj^ indicates that this truncation scheme 
provides reasonable numerical accuracy. The values of the kinematic moment of inertia J^^^ , 
charge quadrupole moment Qo, mass hexadecapole moment Q40, binding energies, separa- 
tion energies S2n{Z,N), and 52n{Z,N) obtained with truncation of the basis at Np — 14 
and Nb = 16 differ from the values obtained with Np = 18 and Nb — 18 by less than 
0.75%, 0.9%, 3.4%, 0.1%, 40 keV, and 40 keV, respectively. The convergence in energy of 
our calculations is similar to that reported in non-relativistic calculations of Ref. [53] based 
on the Gogny force. 

III. ROTATIONAL RESPONSE 
A. The 2^^No ground band 

The observed moments of inertia of the ground band in ^^^No [54-56] are compared 
with the calculated values in Fig. 1. The CRMF calculations without pairing, based on 
the NLl parameterization, marked as CRMF(NLl) in Fig. la, provide an almost constant 
kinematic moment of inertia J^^^ !v 87 MeV"^ up to ~ 0.26 MeV and a dynamic moment 
of inertia which slightly increases with rotational frequency. A band crossing with another 
configuration takes place at ~ 0.26 MeV. These calculations provide a reference point for 
how much the moments of inertia decrease due to pairing. It is interesting to note that the 
moments of inertia in the calculations without pairing are only one half of the rigid body 
value. This unexpected result will be discussed in detail in forthcoming article [45] . 

The CRHB calculations without particle number projection (scaling factor / = 1), 
marked as CRHB (NLI+DIS) in Fig. lb, agree very well with experiment up to = 0.18 
MeV. At higher frequency, experiment and theory diverge. With approximate particle num- 
ber projection using the LN method, marked as NLI-I-DIS-I-LN in Fig. Ic, the theory un- 
derestimates the experimental kinematic and dynamic moments of inertia by ~ 25%. This 
result is in contrast with the good agreement obtained by the same method for supcrde- 
formed bands with A ~ 190 [29,30] and for nuclei in the rare-earth [46] and A ~ 75 [45] 
regions. 
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Different parametrizations of the RMF Lagrangian give quite similar results for the 
moments of inertia if / = 1. For example, the results of the CRHB+LN calculations based 
on the NL3 parametrization [49] (marked as NL3+D1S+LN in Fig. Ic) provide moments of 
inertia which are only shghtly lower (by ^ 3 MeV"-*^) than the ones obtained with the NLl 
parametrization. The CRHB+LN calculations with the NL-Z [48], NLSH [50] and NL-RAl 
[19] parametrizations of the RMF Lagrangian give results which are quite similar to those 
obtained with NL3 ^ and, thus, they are not shown in Fig. 1. The moments of inertia are 
very similar despite the differences in the single-particle spectra near the Fermi level (see 
Figs. 2 and 3). Hence, the most likely reason for discrepancies between experiment and the 
CRHB+LN calculations lies in an inadequate parametrization of the Gogny force (DIS) in 
the particle-particle channel for this mass region, which gives too strong pair correlations. 
This is not quite imexpected since no experimental data above ^°^Pb have been used when 
the DIS set was fitted. The study of other heavy nuclei around ^^"^No also shows that the 
kinematic moments of inertia obtained in the CRHB+LN calculations with the original DIS 
force (scahng factor / = 1.0) are systematically lower than experimental ones (by ~ 20% in 
even-even ^'^^~^^^Pu nuclei). Different parametrizations of the RMF Lagrangian give similar 
results and the deformations of these nuclei are well described in the calculations (see Sect. 
IV). It is unlikely that other available parametrizations for the Gogny force such as Dl [61] 
and DIP [62] will improve the situation, since they produce even stronger pairing than the 
DIS set in superdeformed bands of the ^4 ~ 190 mass region [29]. 

The results of the calculations for the ^^'^No rotational band obtained in the non- 
relativistic cranked HFB approach based on the Gogny force with DIS set of parameters 
[65] seem to support this interpretation. These calculations, which are performed without 
particle number projection, also come very close to the data. One presumes that the in- 
clusion of particle number projection by means of the LN method will lower the calculated 
kinematic moment of inertia, as has been seen in the rare-earth region [46], leading to a 
similar situation as described above. 

B. Selection of the pairing strength 

Quantitative information on the strength of the pair correlations can be extracted from 
the odd-even mass differences, excitation energies of high-K isomers or the moments of 
inertia. We use the moments of inertia for an adjustment of the strength of the Gogny force 
because they are not too sensitive to the details of the single-particle spectrum (see above). 

Our CRHB+LN calculations indicate that in the A ~ 250 mass region the strength of 
pairing correlations should be reduced in order to reproduce the observed moments of inertia. 
The scaling factor / of the Gogny DIS force (see Eq. (10)) has been chosen to reproduce 
the experimental kinematic moment of inertia of ^^^No at rotational frequency Q,x — 0.15 



-"^Thc CRMF calculations (without pairing) for the ground state band in ^^^No also show a weak 
dependence of the moments of inertia on the RMF parametrization, and, thus on details of the 
single-particle structure. A similar situation has been encountered earlier in the ~ 60 [57] and 
A ~ 150 [58-60] mass regions of superdeformation. 
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MeV. The values found for the various parametrizations of the RMF Lagrangian are given 
in Table I. These scaling factors, which are nearly the same, are used in all subsequent 
calculations, unless otherwise specified. The scaled CRHB+LN calculations reproduce the 
amplitude and the f2a;-dependence of the dynamic and the kinematic moments of inertia in 
^^^No (see Fig. 4b,d) rather well. With NL3, experiment and theory agree very well, while 
with NLl some discrepancy develops above Qx = 0.2 MeV. Our choice of the scaling factor 
/ leads also to a reasonable description of the odd-even mass differences in the CRHB+LN 
calculations (see columns 5 and 6 in Table II). 

The need for attenuation of the DIS force within the framework of the CRHB+LN theory 
is not surprising since its pairing properties were adjusted by fitting only the odd-even mass 
differences of the Sn isotopes. Thus the quality of the description of pairing may deteriorate 
far from this mass region. Indeed, the moments of inertia of nuclei in mass regions closer to 
the Sn region, such as the rare-earth region [46], the superdeformed A ~ 190 mass region 
[29,30] and neutron-deficient A ~ 75 region [45], are described well by means of CRHB+LN 
calculations using the original DIS force. 

CRHB calculations (without LN) with original scaling factor / = 1.0 provide a reasonable 
description of both moments of inertia before band crossing (see Fig. lb) and odd-even 
mass differences (see Table II). This approach will be applied to the calculations of the 
quasiparticle spectra in odd nuclei (see Sects. VIE and VI D), for which the CRHB+LN 
calculations are numerically less stable. However, it is not justified for the calculations at 
large rotational frequencies because an unphysical pairing collapse takes place above the 
crossing between the ground and 5'-bands. 

C. High-spin behavior 

Alignment and backbending features of the rotational bands in the actinide region have 
been discussed in a number of publications; see Refs. [66,67] and references quoted therein. 
In this mass region, two high-j shells (^13/2 for protons and J15/2 for neutrons) come close 
to the Fermi surface and the angular momentum of quasiparticles in either orbital can align 
with the axis of rotation. The CRHB+LN calculations with the NL3 parametrization show 
that the alignment of the proton ^13/2 pair (7r[633]7/2) and neutron ji5/2 pair (i/[734]9/2) (see 
Fig. 5) takes place simultaneously in ^^^No at Qx ~ 0.32 MeV (see Fig. 6). The total angular 
momentum gain at the band crossing is ~ 17^, with proton and neutron contributions of 
~ 7h and ~ 10^, respectively. The alignment leads to a decrease of the mass quadrupole 
moment Qq, to a sign change of the mass hexadecapole moment Q^q and to an appreciable 
increase of 7-deformation (see Fig. 7a,b,c). A similar situation holds also in ^^^No, but the 
total angular momentum gain at the band crossing is smaller (r^ llh). 

The simultaneous alignment of proton and neutron pairs occurs also in calculations with 
NLl. The crossing frequency is shifted down by ~ 0.01 MeV (see caption of Fig. 4), the 
angular momentum gain at the band crossing is slightly different, and the high-j particles 
align in ^^^No more gradually. 

Our results differ from those of the cranked HFB calculations based on the Gogny force 
[65], which indicate in ^^^No upbending at / ~ 30^ and backbending at / ~ 38^. These 
calculations are performed without particle number projection, which result in a collapse of 
neutron pairing correlations at relatively low spin, / ~ 20^. In our calculations, the pairing 
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energies decrease with increasing rotational frequency due to the CorioUs anti-pairing effect, 
but there is no collapse of pairing (see Fig. 7d). The experimental data do not extend up 
to predicted backbending and thus do not discriminate between these calculations. 

D. 252no versus ^^^Nq 

The kinematic and dynamic moments of inertia of the band in ^^^No show similar trends 
in rotational frequency as the experimental data (see Fig. 4a, c). In experiment, the moments 
of inertia at low rotational frequencies are smaller for = 150 than for N = 152, in contrast 
with the calculations. One possible reason is the fact that the CRHB+LN calculations give 
deformed shell gaps at iV = 148 and/or 150 (dependent on parametrization), rather than at 
N — 152 as seen in experiment (see Sects. V and VI for details). Indeed, in the calculations 
with the NL3 set the neutron contribution to the total moment of inertia (see Fig. 8) increases 
at the = 148 shell gap, most likely due to the weakening of the neutron pair correlations 
(see the pairing energies for the Fm isotopes in Fig. 9, which are similar to those for the No 
isotopes). This suggests that if the calculations were to give a shell gap at iV = 152, the 
relative magnitudes of the moments of inertia in 252,254]\^q would be reproduced. 

E. Results for other nuclei and general trends of the moments of inertia as functions 

of the particle number 

Figures 10 and 11 demonstrate that the A^-dependence of the moments of inertia in 
254,256pj^ and 248,25o,252q£ jg rather well described in all parametrizations, whereas some 
problems exist for 242-250Qj^ ^ggg pjg Ha). The absolute values are typically reproduced 
within a few % in the Fm and Cf isotopes, but the discrepancy between experiment and 
calculations becomes somewhat larger for the Cm isotopes. The experimental values of 
the moments of inertia in the N = 152 isotopes are reproduced rather well in all RMF 
parametrizations, with the exception of NLSH, which somewhat underestimates the mo- 
ments of inertia (see Fig. 12). One should note, however, that the maximum value of J^^^ is 
at Z — 96 in the calculations, while available data show the maximum ai Z — 98. 

It is interesting to compare the present CRHB+LN results with those from other models. 
The calculations using the MM method in Ref. [68] agree reasonably well with experimental 
excitation energies of E{2'^) states and thus with the moments of inertia (see Figs. 10, 
11 and 12) in this mass region. This is not surprising considering that these data have 
been used in the fit of the strengths (monopole and isospin-dependent) of the proton and 
neutron pairing correlations. However, despite the use of 4 adjustable parameters for pairing 
and better single-particle spectra obtained in the Woods-Saxon potential, the overall level 
of agreement is comparable with that obtained in the CRHB-I-LN calculations (see Figs. 
10, 11 and 12). For example, a detailed examination indicates that the A'"-dependence is 
not described correctly for 248,250Qf 244,246q^ ^g^g rp^^^g ^ ^lei. [68] and Fig. 11 in 
the present manuscript). A similar problem exists for the Z-dependence of the moments 
of inertia in ^^^No and 250Qf nuclei (see Fig. 12). This suggests that the description of 
the single-particle states within the Woods-Saxon potential, with the "universal" set of 
parameters [69] , is still not completely correct in this mass region, despite the fact that the 
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systematics of the experimental data on both spherical and deformed odd-mass nuclei were 
simultaneously taken into account in the fit of these parameters. 

While there are several calculations of the moments of inertia by means of the MM 
method in the actinide and trans- actinide regions (see Ref. [68] and references therein), 
httle has been done in microscopic approaches so far. The rotational bands in ^^^'^^^No have 
been studied in the cranked HFB approaches based on the Skyrme [70,71] and Gogny forces 
[65,72] (see discussion above). The relative magnitude of the moments of inertia of these 
two nuclei is not reproduced in either approach (see, for example. Fig. 6 in Ref. [70]). In the 
calculations with Skyrme forces, this is most likely due to the fact that the deformed shell 
gap appears a,t N — 150 [28], instead of the experimentally observed value oi N — 152. 

The systematic calculations for the Fm isotopes and N = 152 isotones permit the fol- 
lowing general observations. Different RMF parametrizations give similar A^- and Z- depen- 
dencies of the moment of inertia (see Figs. 10 and 12). In the Fm isotopes, the correlation 
between the calculated quadrupole deformations /32 (see Fig. 14) and moments of inertia 
J^^^ (see Fig. 10) is clearly seen with nearly constant values of J^^^ and /32 at A?" = 138 — 160, 
followed by a pronounced drop of both for iV > 160. 

The situation is more complicated in the = 152 isotones, where the change of J^^^ as 
a function of proton number does not correlate with (32- While is almost constant for 
90 < Z < 110 (see Fig. 13), J^^-* shows a maximum at Z = 96 and a minimum at Z = 108 
(see Fig. 12). 

By comparing the different RMF parametrizations, one can see that the sets which 
produce a smaller quadrupole deformation also produce a smaller moment of inertia. NLl 
and NL-Z give very similar J^^\ The same holds for NL3 and NL-RAl, whereas NLSH 
provides smaller values of J^^\ It can also be seen that differences in the underlying single- 
particle spectrum (see Sect. VI) do not lead to significant modifications of the moments of 
inertia. This suggests that many orbitals contribute to the angular momentum. 

F. Summary for Section III 

In general, the moments of inertia for the heaviest nuclei are well described by the 
CRHB-I-LN theory. However, it was necessary to reduce the strength of the DIS Gogny 
force in the pairing channel by ~ 12%, whereas in lighter nuclei with A f=:i 70 — 190, the 
moments of inertia are well reproduced with a full strength DIS force. This points to a 
somewhat different A-dependence of the pairing strength than given by the Gogny force. 
The trends around A ~ 250 with respect to neutron or proton numbers are reasonably 
well reproduced. The calculations with and without pairing indicate very weak dependence 
(within 5%) of the moments of inertia on the parametrization of the RMF Lagrangian. 
The moments of inertia in this mass region are highly collective. Since many single-particle 
orbitals contribute, they are insensitive to fine details of the single-particle states. On the 
other hand, deformed shell gaps affect moments of inertia to some extend leading to a larger 
values. 
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IV. DEFORMATIONS 



A. Comparison with experiment 



Direct experimental information on the deformations of nuclei from Coulomb excita- 
tion and lifetime measurements is quite limited [73]. An alternative method is to derive a 
quadrupole moment from the 2"*" — > 0^ transition energy by employing the relation given by 
Grodzins [74] or by later refinements [73,75]. The prescription of Ref. [75] has an accuracy of 
about 10%. In Figs. 13, 14 and 15, the results of the CRHB+LN calculations are compared 
with deformations extracted by this method. Prom the calculated and experimental charge 
quadrupole moments Q, we derive the deformation parameters P2 by the relation 



The simple linear expression is used to maintain consistency with earlier papers [73]. It is 
sufficient for comparison between calculations and experiment because the same relation is 
used. Including higher powers of /32, e. g. as in Ref. [76], yields values of P2 that are f=:i 10% 
lower. 

Figures 13, 14 and 15 demonstrate that the values of (32 obtained from the 2+ —>■ 0"*" 
transition energies with the prescription of Ref. [75] are consistent with those from Coulomb 
excitation measurements. The CRHB+LN calculations with the NL3, NL-RAl, NLl, and 
NL-Z paramctrizations agree rather well with these values. Considering the uncertainties 
on the extracted values of (32 and the limited experimental data, it is difficult to give any 
preference for a particular set. Only NLSH seems to systematically underestimate (32- 



Figures 13 and 14 illustrate the general trends of deformation as functions of proton and 
neutron numbers for the = 152 isotones and for the Fm {Z = 100) isotopes. In the Fm 
chain, (32 increases gradually from = 138 up to f» 150 for all paramctrizations except 
NLSH, which gives a slight decrease around N = 140. For N > 150, there is a gradual 
decrease of the (32 values, which becomes more rapid above ^ 160. These trends are more 
pronounced with NL3, NLSH and NL-RAl, which have been fitted to the data on neutron- 
rich nuclei as well. The variations arc more gradual in the NLl parametrization, which was 
obtained by fit to data from the beta-stability valley. The mass hexadecapole moments Q40 
are similar for all paramctrizations and decrease with increasing neutron number. While the 
changes of the slope of (32 as a function of neutron number (see Fig. 14a) correlate with the 
shell gaps at A^ = 148, 150 and A^ = 160, 162 (Sect. V), no such correlations are seen for the 
Q40 values (sec Fig. 14b). 

The (32 and Q40 values change more gradually with proton number Z in the A^ = 152 
chain (see Fig. 13). The (32 values are almost constant as a function of proton number. 
The calculated mass hexadecapole moments Q40 show a sinusoid-like curve as a function of 
proton number, with a maximum at Z 94 and a minimum ai Z ^ 106. 

We expect that the trends are similar in the chains adjacent to the Fm and A^ = 152 
chains, which is corroborated by the less systematic calculations for the Cm, Cf and No 




(14) 



B. General trends 
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isotopes (see Fig. 15). The equilibrium deformations are very similar for NL3 and NL- 
RAl on the one hand, and for NLl and NL-Z on the other hand. For this reason, the 
results obtained with NL-RAl and NL-Z are omitted in Figs. 13, 14 and 15. The calculated 
/?2 increases as the parametrization changes from NL-SH to NL3 to NLl. This trend has 
previously been seen in the A ~ 60, 150, and 190 regions of superdeformation [30,57,59]. 

C. Summary for Section IV 

In summary, the CRHB+LN theory with the NL3, NLl, NL-RAl and NL-Z parameter 
sets satisfactorily reproduces the magnitude of the /?2 deformation of the heaviest nuclei, 
where they have been measured, whereas NLSH systematically underestimates it. 

V. SHELL STRUCTURE 

The stability of the superheavy elements is due to a region of low level density in the 

single-particle spectrum. For deformed nuclei all single-particle states are two-fold degener- 
ate, and, thus, the region of low level density generally correlates with the 'shell gap'. The 
situation is more complicated in spherical nuclei, where the shell correction energy depends 
not only on the size of the shell gap, but also on the degeneracy of single-particle states in the 
vicinity of the Fermi level [77,78] . It is a concern of this paper to study how well the different 
RMF parametrizations reproduce the shell gaps in the heavy deformed elements. A simple 
intuitive measure for the level density at the Fermi surface is the distance Esp-gap between 
the last occupied and the first empty levels. Another way is to consider the two-neutron 
S2niZ,N) and two-proton S2p{Z,N) separation energies 

S2n{Z, N) = B{Z, N) - B{Z, N-2) 

S2j>{Z, N) - B{Z, N) - B{Z - 2, AT), (15) 

where B{N,Z) is the binding energy. The separation energies show a sudden drop at the shell 
gaps, if they are large. If the variations of the level density are less pronounced, the quantity 
52n{Z.,N) related to the derivative of the separation energy is a more sensitive indicator of 
the localization of the shell gaps. For the neutrons (and similarly for the protons) it is 
defined as 

82n{Z, N) = S2n{Z, N) - S2n{Z, N + 2) = 

= -B{Z, N-2) + 2B{Z, N) - B{Z, N + 2). (16) 

We show in Appendix A that variations (but not their absolute values) of 62n{Z,N) and 
Esp-gap agree rather well. 

In this section, we study the shell structure along both the Fm (Z — 100) and the 
= 152 chains. 
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A. Shell structure along the Z = 100 line 



The results for two- neutron separation energies S2n{Z, N) in the Fm (Z = 100) chain 
with different RMF parametrizations are shown in Fig. 16. There is a systematic difference 
between the NLl, NL-Z sets and the NL3, NLSH, NL-RAl sets. The former underestimate 
two-neutron separation energies, thus reveahng their weakness in the description of isotopic 
trends. NL-Z is somewhat better as compared with NLl. Considering that these two sets 
are fitted to the same data, this resuh together with the one shown in Fig. 2 of Ref. [52] 
possibly indicate the importance of a more microscopic treatment of the center-of-mass 
correction for the reproduction of isotopic trends. On the other hand, NL3, NLSH and 
NL-RAl better reproduce the experimental S2n{Z. N) values (see Fig. 16) indicating that 
the isovcctor component of the interaction has been treated more carefully in these sets. 

The latter sets reasonably describe S2n{Z,N) for 144 < TV < 148 and 154 < iV < 159, 
but underestimate S2n{Z,N) for N — 150,152. The shoulder in experimental S2n{Z,N) 
values at N — 152 reveals a deformed shell gap [79], which is better seen in the plot of 
S2n{Z,N) (Fig. 17a). The size of this gap depends sensitively on proton number. As seen 
in Fig. 17, there are discrepancies with experiment: NL3 and NL-RAl (NLl and NL-Z) 
produce a gap at iV = 148 (iV = 148, 150) instead of at iV = 152 and NLSH does not show 
a clear gap (see also Figs. 3 and 19). The analysis of the neutron quasi-particle spectra in 
249,25iQf ^ggg Sect. VI D) also indicates that the calculated shell gaps do not correspond to 
the experimental ones. 

Earlier calculations predicted the presence of a deformed neutron shell gap at = 162 
in superheavy nuclei (see, for example, Refs. [14,28] and references quoted therein). The 
presence of this gap in nuclei with Z fa 108 was confirmed by recent experimental information 
[10,12]. The appearance of this gap and its size strongly depend on the parametrization of 
the specific theory and on the proton number. For example, in the Skyrme Hartree-Fock 
calculations [28], this gap is pronounced in the Skl3 parametrization, where it is seen over 
the proton range of Z = 98 — 116, but is absent in the SkP parametrization. This gap is 
clearly seen in RMF calculations with the NL3 and PL-40 parametrizations [28] but only 
at proton numbers around Z — 106. The present CRHB-I-LN calculations in the Fm chain 
indicate a gap at A^ = 162 for NL3 and at A^ = 160 for NLl, NL-Z and NL-RAl (see Fig. 
17). However, the small value of 52n{Z, N) 0.8 MeV suggests a small gap. 

B. Shell structure along the N = 152 line 

To judge the reliability of predictions of superheavy nuclei, it is critical to see how differ- 
ent RMF parametrizations are able to describe the experimental shell gaps as a function of 
proton number. The calculations for the A^ = 152 isotones are compared with experimental 
data in Figs. 20 and 21. As for the two-neutron separation energies, one can see that the 
two-proton separation energies S2p{Z,N) arc best described by NL3, NLSH, NL-RAl. In 
contrast, the S2p{Z, N) values are overestimated by NLl and NL-Z, which were obtained 
by fit to stable nuclei. The S2p{Z, N) plots do not show clearly where the proton deformed 
gaps are located, which becomes visible in the 52p{Z, N) plots. The experimental data show 
a shell gap aX Z = 100. Only NLSH describes the position of this gap and the 52p{Z, N) 
values agree very well. However, the analysis of the quasi-particle spectra in Sect. VI G 
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reveals that this gap hes between the wrong bunches of single-particle states. Calculations 
with NLSH also indicate a gap at Z = 108, which has not been observed so far. NL-RAl 
does not show any deformed gap for 92 < Z < 108 (see Fig. 21). NL3, NLl and NL-Z 
give a shell gap at Z — 104, in contradiction with experiment. The analysis of the proton 
quasi-particle spectra in ^^^Bk (see Sect. VI E) leads to the same conclusion. 

Many effective interactions not specifically fitted to the actinide region encounter similar 
problems. For example, most of the Skyrme forces fail to reproduce the deformed Z = 100 
shell gap in the N = 152 isotones (see Fig. 5 in Ref. [28]). SkI4 is the only force which 
shows this gap. The SkI3, Skll and Sly6 forces show a Z — 104 sheU gap, while the SkM* 
and SkP forces do not show any gap sit Z — 100 — 104. 

C. Pairing along the Z = 100 and N = 152 lines 

Figure 9 shows the pairing energies Epairing (sec Eq. (11)) obtained with NLl and NL3 
as a function of neutron number along the Z = 100 line and as a function of proton number 
along the N = 152 line. The general trend as a function of nucleon number is the same for 
both sets. 

Let us first discuss the Z = 100 line. In both RMF parametrizations, neutron pairing 
correlations arc weakest at ^ 148, reflecting the presence of a shell gap at = 148 in 
the NL3 paramctrization and somewhat smaller gaps at = 148 and A^ = 150 in the NLl 
parametrization (see Fig. 3). Going away from these shell gaps, the neutron pairing energies 
increase in absolute value, reflecting the increasing level density (Fig. 19). The neutron 
pairing is weakened at A^ ~ 160 due to the presence of smaller shell gaps at N — 160 (NLl, 
see Sect. VA) and at A^ = 162 (NL3, see Fig. 19). The weakening of the neutron pairing at 
A^ ^ 148 and A^ ~ 160 is more pronounced in the NL3 parametrization as compared with 
NLl, reflecting larger shell gaps (see Fig. 3). 

The proton pairing shows the same trend as the neutron pairing, but with much smoother 
changes in neutron number. In both parametrizations, proton pairing is smaller and stays 
relatively constant at A^ = 138 — 150, reflecting low and nearly constant level density below 
the Fermi level (see Figs. 2 and 22). For A^ > 150, the deformation modifications induced 
by changes in the neutron number increase the proton level density near the Fermi level (see 
Fig. 22), enhancing the proton pairing. 

The pairing energies on the A^ = 152 fine exhibit the same features as discussed above. 
Both in the neutron and proton subsystems they reflect the presence of the Z — 104 and 
Z = 96 shell gaps (see Figs. 2 and 9). 

D. Summary for Section V 

The experimental 52n{Z,N) quantity shows a distinct deformed shell gap at A" = 152, 
which is quite pronounced for No and Fm nuclei and less so for Cm and Cf (see Fig. 18). 
For the Fm isotopes, the CRHB+LN calculations predict a deformed shell gap(s) at A^ = 
148 (NL3 and NL-RAl) or at A^ = 148,150 (NLl and NL-Z); see Fig. 17. The NLSH 
parameterization does not give a clear gap. The experimental data for the A^ = 152 isotones 
show a shell gap at Z = 100; NLl, NL-Z and NL3 give a gap at Z ^ 104 and NL-RAl no 
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gap at all. Only NLSH gives a gap at Z = 100; however, it lies between the wrong bunches 
of single-particle states (see Sect. VI). This demonstrates the fact that the usual analysis 
of shell structure, in terms of only S2n{Z, N) and/or d2n{Z, N), maybe insufficient to judge 
the quality of the parametrization. 

VI. QUASI-PARTICLE STATES 

The investigation of the single-particle states in the A ~ 250 deformed mass region can 
shed additional light on the reliability of the RMF predictions of the energies of spherical 
subshells responsible for 'magic' numbers in superheavy nuclei. This is because several 
deformed single-particle states experimentally observed in odd nuclei of this mass region 
(see Table III) originate from these subshells. Considerable deviations between experiment 
and theory for a specific deformed state will indicate that the position of the spherical 
subshell from which this state originates is poorly described. 

In the past, the RMF studies of single-particle spectra have been mostly performed in 
spherical or near-spherical nuclei (see, for example, Ref. [24] and references quoted therein), 
where a number of restrictions, such as the neglect of the currents or of the breaking of 
the Kramer's degeneracy, have been imposed in order to simplify the task. In addition, 
the BCS approximation was used. Moreover, a direct comparison between experimental 
and theoretical single-particle states in spherical nuclei should include the particle- vibration 
coupling, which can modify the single-particle energies considerably [27]. The modification 
of the quasiparticle states by particle- vibration coupling is weaker in deformed nuclei [84-86] . 

Not much is known about the accuracy of the description of the quasiparticle states 
in deformed nuclei within the framework of the RMF theory. In most cases the analysis 
of odd nuclei was based on the single-particle spectra calculated in neighboring even-even 
nuclei. To our best knowledge, a direct comparison between experimental and theoretical 
quasi-particle spectra has not been published. 

The most detailed attempt to analyze the single-particle properties at finite deformation 
and their connections with those at spherical shape has been done within the cranked RMF 
theory in Ref. [59] for superdeformed (SD) bands in the A ~ 150 mass region. Although 
a direct comparison between experimental and calculated single-particle energies was not 
possible, some general conclusions were drawn based on a systematic analysis of experimental 
data and the expected response of specific single-particle states to a change of deformation. 
It was found that the RMF theory provides a reasonable description of the single-particle 
states in the vicinity of the SD shell gaps. However, some deviations between experiment and 
theory were detected, which could reach pa 1 MeV for some states. For example, the relative 
positions of the i/[651]l/2 and i/[642]5/2 states from the u2gg/2 and i^lii3/2 spherical subshells 
are not reproduced. This problem exists in the NLl, NL3 and NLSH parametrizations of 
the RMF Lagrangian. 

A. Computational details 

In the present manuscript, we address for the first time the question of a fully self- 
consistent description of quasiparticle states in the framework of the RMF theory. A proper 
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description of odd nuclei implies the loss of the time-reversal symmetry of the mean-field, 
which is broken by the unpaired nucleon. The BCS approximation has to be replaced by 
the Hartree-(Fock-)Bogoliubov method, with time-odd mean fields taken into account. The 
breaking of time-reversal symmetry leads to the loss of the double degeneracy (Kramer's 
degeneracy) of the quasiparticle states. This requires the use of the signature or simplex basis 
in numerical calculations, thus doubling the computing task. Furthermore, the breaking of 
the time-reversal symmetry leads to nucleonic currents, which cause the nuclear magnetism 
discussed in Sect. II A. 

The CRHB(-I-LN) theory takes all these effects into account. Thus, the effects of blocking 
due to odd particle are included in a fully self-consistent way. The CRHB computer code is 
set up in a signature basis and in three-dimensional Cartesian coordinates. The latter allows 
one to study also the 7-deformation. In order to specify the detailed structure of blocked 
orbitals, the existing CRHB code [30] has been extended to describe odd and odd-odd 
nuclei. The blocking procedure is implemented according to Refs. [26,87,88]. The blocked 
orbital can be specified either by its dominant main oscillator quantum number N or by 
the dominant VL quantum number (f2 is the projection of the total angular momentum on 
the symmetry axis) of the wave function, or by combination of both. In addition, it can be 
specified by the particle or hole nature of the blocked orbital. Note that f2 is not a conserved 
quantum number in the CRHB code. As a consequence, convergence problems, emerging 
from the interaction of the blocked orbital with others, appear somewhat more frequently 
than in a computer code restricted to axial symmetry. Convergence problems appear more 
frequently when approximate particle number projection by means of the Lipkin-Nogami 
method is imposed, which is most likely due to additional non-linearities. 

As illustrated in Fig. 23, the quasi-particlc spectra calculated within the CRHB-I-LN 
(with LN) framework and using the scaled DIS Gogny force given in Table I are very similar 
to those obtained by means of the CRHB (without LN) with the original DIS force. The 
difference in the energies of the quasi-particle states is typically less than 150 keV and 
the level ordering is the same. Thus in order to avoid the convergence problems in the 
calculations with LN, all other calculations of quasi-particle states were performed in the 
CRHB framework with the original DIS force. 

B. Particle- vibration coupling and other effects 

Figures 23, 24 and 25 show that the calculated quasiparticle spectra are less dense than 
in experiment. The average level density of the single-particle states is related to the ef- 
fective mass (Lorentz mass in the notation of Ref. [90] for the case of RMF theory) of 
the nucleons at the Fermi surface ni*{k},^)/m. The RMF theory gives a low effective mass 
m*{kF)/m 0.66 [17]. The experimental density of the quasiparticle levels corresponds to 
an effective mass m*{kF)/fn close to one. This discrepancy appears also for non-rclativistic 
mean-field models [23]. It has been demonstrated for spherical nuclei that the particle- 
vibration coupling brings the average level density in closer agreement with experiment [27], 
which means m*{kF)/m closer to one. In a similar way, the particle- vibration coupling leads 
to a compression of the quasi-particle spectra in deformed nuclei [86]. The surface vibra- 
tions are less collective in deformed nuclei than in spherical ones because they are more 
fragmented [84,91]. As a consequence, the corrections to the energies of quasiparticle states 
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in odd nuclei due to particle-vibration coupling are less state-dependent in deformed nuclei. 
Hence the comparison between experimental and mean field single-particle states is less 
ambiguous in deformed nuclei as compared with spherical ones [27,91], at least at low exci- 
tation energies, where vibrational admixtures in the wave functions are small. Calculations 
within the quasiparticle-phonon model [92,93] indicate that in the A ~ 250 mass region the 
lowest states have mainly quasi-particle nature and the corrections to their energies due to 
particle-vibration coupling are typically at the level of 150 kcV or less. The states above 
w 700 keV contain very large vibrational admixtures [86] and thus experimental states above 
this energy should not be compared with the pure quasiparticle states obtained in the CRHB 
calculations. 

Since particle-vibration coupling is not included, it is important to estimate how large is 
the discrepancy between calculated and experimental quasiparticle energies due to the low 
effective mass m*{kF)/m of the RMF theory. Assuming for an estimate that the effective 
mass just stretches the energy scale, the difference between the energies of quasiparticle 
states obtained in the calculations with ^^^^ and = 1 is 

m m 

AE, = E„i'^)(l-'^). (17) 
m m 

which remains below ~ 200 keV as long as the calculated state is located in the energy 
window of 700 keV with respect to the Fermi surface, whereas it grows for higher excitation 

energies. 

We are mostly interested in how well the positions of spherical subshells are described 
in the RMF calculations. One may reduce the error by comparing the experimental and 
calculated averages of two (or more) deformed single-particle states emerging from the same 
spherical subshell. The average of deformed states has the same energy when extrapolated 
to spherical shape as each of these states. The advantage is that the average of the states 
has smaller excitation energy than at least one of these states. As a result, the energy of the 
spherical subshell can be estimated more precisely (within 200 keV if \ [Eqp^i+Eqp^2) /'^\ < 700 
keV). Such an approach is especially useful when one of the states has particle {E'^p > 0) 
and the other hole {E^^ < 0) nature, since their average energy can be well within the 700 
keV energy window even if the excitation energy of each is far outside this window. 

An additional source of uncertainty is the Coriolis interaction between the quasi-particle 
states, which is neglected. It is relatively modest, affecting the energies of quasi-particle 
states by at most 100 — 200 keV [63,94]. As a whole, the uncertainty of our estimate for the 
spherical subshell energies in Sects. VI E and VI D is around 300 keV. 

C. General observations 

Figures 2 and 3, where the single-particle states in ^^^No are shown for different 
parametrizations, reveal important trends. The energies of some single-particle states de- 
pend strongly on the parametrization. For example, for 7r[521]3/2 and 7r[521]l/2 ^, the 



^These states arc of special interest since they originate from the spin-orbit partner spherical 
subshells i^lj-jji and 'Klj^jii which define the size and the position of the magic spherical proton 
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single-particle energies calculated with NLl and NLSH differ by ?a 2 MeV. The small differ- 
ences in the self-consistent deformations cannot explain the differences in the single-particle 
energies. 

Another observation is that the relative energies of the different Q states, which emerge 
from the same spherical subshell, almost do not depend on the parametrization. For protons, 

these are, for example, the 7r[642]5/2, 7r[633]7/2, 7r[624]9/2 states from the 7rlii3/2 spherical 
subshell and the vr [5 14] 7/2, tt [505] 9/2 states from the vr 1/19/2 subshell. This is expected 
because the splitting of different fl states from the same spherical subshell is a consequence 
of the deformation of the mean field, which is not very sensitive to the parametrization (see 
Sect. IV). On the other hand, the single-particle energies of the 712/7/2 and 712/5/2 spherical 
subshells (but not their splitting), as well as the deformed states emerging from them, change 
considerably when going from the NLl to the NLSH parametrization (see Fig. 20 in Ref. 
[59] and Fig. 2 in the present manuscript for their deformed counterparts). This leads to a 
Z — 104 deformed shell gap in the NLl and NL-Z parametrizations and a, Z — 100 gap in 
the NLSH parametrization. 

The situation is analogous for the neutron states (see Fig. 3). For example, the rela- 
tive energies of the z/[624]7/2, z/[615]9/2 and z/[606]ll/2 states, emerging from the z/liii/2 
spherical subshell, almost do not depend on the RMF parametrization. The same is true 
for relative energies of the states i/[743]7/2, z/[734]9/2 and i/[725]ll/2, emerging from the 
i^lji5/2 subshell, and the states i/[761]l/2* and z/[752]3/2* from the z^lji3/2 subshell. It is 
interesting that the relative energies of the states z/[631]l/2 {v2g'ji2), z/[622]5/2 {i>2gg/2), 
i/[622]3/2 (z/2^7/2), z/[620]l/2 (2/3^5/2), z/[613]7/2 (2/2^9/2) and z/[613]5/2 (2/2^7/2), originat- 
ing from the different spherical subshells shown in parenthesis after the Nilsson labels, are 
almost independent on the RMF parametrization. This indicates that the relative energies 
of the v2gg/2i i^2g-j/2 and 1^3^5/2 spherical subshells only marginally depend on the RMF 
parametrization, which is clearly seen in the single-particle spectra of spherical nuclei (see 
Fig. 20 in Ref. [59] and Fig. 1 in Ref. [17] for spectra of ^°^Pb and Figs. 4, 9 and 15 in 
Ref. [17] for spectra of superheavy nuclei). On the other hand, the increase of the sepa- 
ration between the z^[761]l/2* and z/[750]l/2* states when going on from the NLl to the 
NLSH parametrizations shows that the separation between the z^lji3/2 and v2hii/2 spherical 
subshells increases. 

In a previous study in the A ~ 150 mass region of superdeformation [59], we also con- 
cluded that the dependence of the single-particle energies of deformed states on the RMF 
parametrizations reflect their energy displacement at spherical shape. Since different RMF 
sets give similar spin-orbit sphttings (see Fig. 2 in Ref. [17] and Fig. 2 for the splitting and 
position of the tt [521] 3/2 and 7r [521] 1/2 states), the dominant modification is a shift of the 
position of the /-shells. 

There are similarities between the single-particle spectra (see Figs. 2 and 3) obtained 
with the NLl and NL-Z parametrizations, on the one hand, and NL3 and NL-RAl, on the 



shell gap in superheavy nuclei [17]. Their splitting, defined primarily by the spin-orbit splitting 
and their response to deformation, almost does not depend on the RMF parametrization. In part, 
this is due to the fact that their interaction with other orbitals is rather weak (see, for example, 
the Nilsson diagram (Fig. 4) in Ref. [79]). 
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other hand. Thus detailed study of quasi-particle spectra will be performed with only the 
NLl and NL3 parametrizations. They can be considered as representative examples of the 
two groups of RMF parametrizations discussed in Sect. II B. 

D. Odd-neutron nuclei 249,251 Qf 

The quasiparticle spectra of these two nuclei with neutron numbers N — 151, 153 are 
presented in Figs. 23 and 24. The ground state configuration of 249Qf jg correctly reproduced 
in both parametrizations. Only NLl gives the correct ground state i/[620]l/2 in ^^^Cf, 
whereas NL3 gives the z/[615]9/2. 

The z/[622]5/2 and z/[613]7/2 (and z/[604]9/2 in ^s^Cf) states emerge from the u2gg/2 
spherical subshell. The i/[622]5/2 energy is reproduced within 300 keV in both nuclei by 
both parametrizations. However, in both nuclei the excitation energy of the z/[613]7/2 state 
is overestimated by ^ 0.55 MeV and hj ^ 1.0 MeV in the NLl and NL3 parametrizations, 
respectively. The comparison of the average energies of the experimental and calculated 
i/[613]7/2 and i/[622]5/2 states suggests that the energy of the i'2gQ/2 subshell has to be 
decreased by f« 0.15 MeV and by ^ 0.6 MeV in the NLl and NL3 parametrizations, respec- 
tively. 

The relative position of the z/[613]7/2 and z/[615]9/2 states, emerging from the i'2gg/2 
and v\i\\i2 spherical subshells, is not reproduced in both parametrizations. In addition, the 
excitation energy of the z/[624]7/2 state from the t^liii/2 subshell is overestimated by > 1 
MeV in 249q£ ^ggg pjg 23). An analysis similar to the one given above suggests that the 
energy of the z/Hii/2 spherical subshell has to be increased by 0.3 MeV and k,\ MeV 
in the NLl and NL3 parametrizations, respectively, in order to bring the calculations in 
agreement with experiment. 

In 25iCf, the relative positions of the i/[622]3/2 {v^gy^) and z/[620]l/2 (z/34/2) states 
(see Fig. 24) differ in experiment and calculations (with NL3). In addition, NL3 fails to 
reproduce the i/[620]l/2 energy. The states i/[622]3/2 and i/[620]l/2 are almost degenerate 
in energy in the NLl parametrization, while in experiment the excitation energy difference 
between the z/[622]3/2 and z/[620]l/2 states is 178 keV. An increase (decrease) of the energy 
of the v^g-jj^ spherical subshell by ~ 0.15 MeV in NLl (NL3) and a decrease of the energy 
of the v'id^i2 spherical shell by 0.5 MeV in the NL3 parametrization would remove these 
discrepancies. 

In addition, the z/[761]l/2* ^ state has been observed in ^^-"^Cf. Its quasiparticle energy 
is well described (within 150-200 keV) in both parametrizations, suggesting that the energy 
of the i^lji3/2 spherical subshell is correctly accounted for in both parametrizations. 

The comparison of the average energies of the experimental and calculated z/ [734] 9/2 and 
i/[725]ll/2 states in ^^^Cf (see Fig. 24), from the v\jxhi2 subshell, suggests that the energy 
of this subshell has to be decreased by ~ 0.7 MeV and by ~ 0.45 MeV in the NLl and NL3 
parametrizations, respectively. However, the z/[734]9/2 state would still remain the ground 
state of Qf after all the modifications discussed above. 



^The use of asterisk at the Nilsson labels is explained in caption of Fig. 2. 
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E. Odd-proton nucleus ^*^Bk 



The quasiparticle spectrum of ^^^Bk [Z = 98, N = 152) is presented in Fig. 25. Three 
states 7r[642]5/2, 7r[633]7/2 and 7r[624]9/2, from the 7rlii3/2 subsheU, have been observed in 
experiment. We select this subshell as a reference, with respect to which the positions of 
other spherical subshclls will be compared, because in the NLl and NL3 paramctrizations the 
TT [633] 7/2 state is close to experiment, and the tt [642] 5/2 and tt [624] 9/2 states are located 
below and above the Fermi level, respectively. The 7r[642]5/2 and 7r[633]7/2 states are 
reasonably well reproduced. The tt [624] 9/2 state is excluded from the direct comparison 
since vibrational admixtures are expected to be large due to high experimental excitation 
energy of this state. 

The energy difference between the 7r[521]3/2 and 7r[521]l/2 states is well reproduced, 
suggesting that the spin-orbit splitting between the 7r/5/2 and 7r/7/2 spherical subshells is 
correctly described. However, their positions with respect to the Fermi level depend on the 
parametrization: for only NLl is the tt [521] 3/2 state lowest in energy, in agreement with 
experiment. The energy difference between the [633]7/2 and [521]3/2 states is small, around 
250 kcV (sec Fig. 25). A decrease (increase) of the energy of the 2/7/2 spherical subshells 
by ~ 0.25 MeV in the NL3 (NLl) paramctrizations would bring the relative positions of 
these states in agreement with experiment. This also implies the same shift of its spin-orbit 
partner 772/5/2- 

The 7r[514]7/2 state, from the 7rl/i9/2 subshell, is too low in energy in both paramctriza- 
tions with respect to the vr [633] 7/2 state. An increase of the energy of the 1/19/2 spherical 
subshell by ~ 0.85 MeV for NL3 and by ~ 0.6 MeV for NLl would bring calculations in 
agreement with experiment. The analyses of odd-proton nuclei around ^o^Pb (see Fig. 7 in 
Ref. [24]) and of shape coexistence in the Pt-Hg-Pb isotopes [95] also point to this deficiency 
in the description of the 1/19/2 spherical subshell energy. 

The calculations underestimate the position of the 7r[530]l/2 state, from the vr2/7/2 sub- 
shell, by ~ 1.5 MeV. If the energy of the l/?.9/2 subshell were increased as discussed above, 
this would push the 7r[530]l/2 state closer to the Fermi level due to the interaction with 
the 7r[541]l/2 state, from the 7rl/i9/2 subshell (see Fig. 4 in Ref. [79]). This would lead to a 
better agreement between calculations and experiment. 

The 7r[400]l/2 state, from the 7r3si/2 subshell, is also reasonably well reproduced (some- 
what better with NL3 than with NLl). However, the relative positions of the 7r[400]l/2 and 
7r[642]5/2 states (the latter from the vr 1^13/2 subshell) suggest that the energy of the 7r3si/2 
subshell has to be increased by fa 0.3 MeV in the NLl parametrization. 

F. Consequences for deformed shell gaps in the A ^ 250 mass region. 

Figure 26 shows how the proton and neutron spectra in ^^^No are modified if the spherical 
subshells are shifted as discussed in Sects. VIE and VI D. Similar corrected spectra are 
obtained with NL3 and NLl, indicating that the shifts are correctly defined. These shifts 
would lead to the deformed shell gaps at = 152 and Z — 100, as seen in experiment 
(Sects. VA and VB). Neutron gaps at A^ = 148,150 (NLl) {N = 148 in NL3) seen in 
uncorrected spectra disappear, while the proton gap at Z — 104 becomes smaller. 
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In addition, the ordering of the neutron and proton single-particle states below and 
above these shell gaps would be more similar to the Woods-Saxon potential (see Fig. 8 in 
Ref. [96]), whose parameters were defined by an overall fit to the single-particle spectra in 
heavy actinide nuclei [97]. 

These examples illustrate that in the region of high level density and small shells gaps, a 
shift of the energies of one or two single-particle states by a modest energy of 0.5 MeV can 
modify the nucleon number of the shell gap by two or four units. A new parametrization 
of the RMF Lagrangian, which implements the shifts discussed in Sects. VIE and VI D 
naturally, is called for. 

G. Estimates for other parametrizations 

The calculations of odd-A nuclei performed with the NLl and NL3 sets indicate that, 
in general, the results are quite similar but somewhat better agreement is obtained in the 
NLl parametrization. Moreover, some conclusions about the accuracy of the description 
of the quasi-particle states in other RMF parametrizations can be drawn with the aid of 
the spectra presented in Figs. 2 and 3. The NL-Z parametrization gives single-particle 
spectra in between those for NLl and NL3 and, thus, a similar accuracy of the description 
of quasi-particle states is expected. 

The agreement with experiment is worse for the NLSH and NL-RAl parametrizations. 
Let us illustrate this by a few examples for NLSH, which deviate most from experiment. 
The energy splitting between [6 13] 7/2 and z/ [6 15] 9/2 increases from Ri 0.5 MeV up to fa 2.5 
MeV when going from NLl to NLSH (see Fig. 3). Thus in order to reproduce the relative 
positions of these states in 249,251 (~;f (^ggg Figs. 23 and 24), the relative distance between the 
spherical i'2gg/2 and z/liii/2 subshells should be corrected by roughly 2 MeV. 

The Fermi level for the odd-proton (Z=97) ^"^^Bk nucleus will be located somewhere in 
the vicinity of the 7r[633]7/2 and 7r[514]7/2 states (see Fig. 3). Thus, the 7r[521]3/2 and 
TT [521] 1/2 states (and the corresponding 712/7/2 and Tilj^ji spherical subshells) should be 
lowered by roughly 1 MeV with respect of the 7r[633]7/2 state (the 7rlii3/2 subshell) in 
order to reproduce the experimental spectra (see Fig. 25). In addition, the position of the 
7r[514]7/2 state (7rl/i9/2 spherical subshell) should be raised by ~ 700 keV with respect of 
the 7r[633]7/2 state (the 7rlii3/2 subshell). 

Thus the empirical shifts required to reproduce experimental quasiparticle energies are 
much larger for NLSH than the ones needed for NLl and NL3 (see Sects. VI D and VIE). 
Only after these shifts will the 7r[633]7/2 and 7r[521]3/2 states be located in the vicinity 
of the proton Fermi level and there will be a gap at Z = 100 between these states and 
TT [521] 1/2 and tt [5 14] 7/2 in agreement with an analysis based on the Woods-Saxon potential 
[96] . Although the NLSH parametrization is the only one parametrization which reproduces 
the Z — 100 gap (see Fig. 21), this gap is created between the wrong bunches of states. 

H. Consequences of nuclear magnetism for quasiparticle states 

The infiuence of nuclear magnetism on the binding energies of one-quasiparticle states 
in 249Qf g^jj^ ^^^Bk is shown in Table IV. In all cases, it provides small additional binding. 
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It is state dependent and lies between —16 and —69 keV, depending weakly on the strength 
of the pair correlations and on particle number projection. This indicates that these heavy 
nuclei are rather robust against polarization effects induced by nuclear magnetism. Thus, if 
nuclear magnetism was neglected, the quasiparticle spectra in this mass region would only 
be marginally modified. The influence of nuclear magnetism on quasiparticle energies is 
larger in lighter systems [45] (see also Ref. [98] for a study of the effects of the time-odd 
mean fields on the quasiparticle energies within the Skyrme Hartree-Fock approach) and in 
two-particle configurations [99]. 

While the neglect of nuclear magnetism seems to be a reasonable approximation for the 
one-quasiparticle energies, it has to be taken into account when the strength of the pairing 
correlations is fitted to experimental odd-even mass differences because it modifies A^^^ by 
^ 10% (see Table II). 

I. Implications for the study of superheavy nuclei 

In the NLl and NL3 parametrizations, the energies of the spherical subshells, from which 
the deformed states in the vicinity of the Fermi level of the A ~ 250 nuclei emerge, are de- 
scribed with an accuracy better than 0.5 MeV for most of the subshells (see Figs. 27 and 
28 where required corrections for single-particle energies are indicated). The discrepancy 
reaches 0.6-1.0 MeV for the vrl/ig/a (NL3, NLl), ulin/2 (NL3), z^lji5/2 (NLl) and u2gg/2 
(NL3) spherical subshells. Considering that the RMF parametrizations were fitted only 
to bulk properties of spherical nuclei, this level of agreement is good. However, the accu- 
racy of the description of single-particle states is unsatisfactory in the NLSH and NL-RAl 
parametrizations (see Sect. VI G). 

The single-particle levels of spherical magic superheavy nuclei are not modified much 
with the empirical shifts of Sects. VIE and VI D (see Figs. 27 and 28 for the calculated 
and corrected single-particle spectra of a ^72^20 nucleus). This conclusion relies on the as- 
sumption that the shifts should be similar in the deformed A ~ 250 mass region and in 
superheavy nuclei. The corrected spectra from the NLl and NL3 calculations are very sim- 
ilar, with minor differences coming from the limited amount of information on quasiparticle 
states used in the analysis. More systematic study of quasiparticle states in deformed nuclei 
are required to determine these corrections more precisely. 

Let us consider the calculations for the nucleus with Z = 120, A?" = 172. The corrected 
single-particle levels still suggest that = 172 and = 184 are candidates for magic 
neutron numbers in superheavy nuclei. The position of the z/4si/2 spherical subshell and the 
spin-orbit splitting of the 3^5/2 and 3d3/2 subshells will decide which of these numbers (or 
both of them) is (arc) actually magic. The corrected proton levels indicate that the Z — 120 
gap is large whereas the Z — llA gap is small. Hence, on the basis of the present investigation 
we predict that Z = 120 is the magic proton number. This conclusion is based on the 
assumption that the NLl and NL3 sets predict the position of the 7rliii/2 and vr3j9i/2,3/2 
subshells within 1 MeV. The positions of 7rlji5/2 and 7T2gQ/2 seems less critical, because they 
are located well above this group of states both in Skyrme and RMF calculations [17]. It 
seems possible to obtain information about the location of the 7rliii/2 subshell, which may 
have been observed through its deformed state (7r[651]l/2*) in superdeformed rotational 
bands of Bi-isotopes [100,101]. An CRHB analysis may provide this critical information. 
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In this context it is important to mention that the RMF parametrizations, NL-SH and 
NL-RAl, which are the only ones to predict Z = 114 as the magic proton number [18,19], 
provide poor descriptions of the single-particle states (see Sect. VI G). 

The Nilsson diagrams given, for example, in Figs. 3 and 4 of Ref. [79], suggest that 
spectroscopic studies of deformed odd nuclei with proton and neutron numbers up to Z 108 
and ~ 164 may lead to observation of the deformed states with Q = 1/2, emerging from 
the vrHii/2 and 7rlji5/2 spherical subshells (located above the Z = 120 shell gap) and from 
i'lkn/2 and either i'2hu/2 or i^lji3/2 subshells (located above the = 184 shell gap). This 
will further constrain microscopic models and effective interactions. 

No information on low-j states, such as 7r3p3/2, 7r3pi/2, 1^^3/2 and i/4si/2, which decide 
whether Z = 120 or Z = 126 and N = 172 or = 184 are magic numbers in micro- 
scopic theories (see Refs. [17,20] and references quoted therein), will come from the study 
of deformed nuclei (see Table III) . 

The measured and calculated energies of the single-particle states at normal deformation 
provide constraints on the spherical shell gaps of superheavy nuclei. In particular, the small 
splitting between the 7r[521]l/2 and 7r[521]3/2 deformed states, emerging from the 7r2/5/2 
and 7r2/7/2 spherical subshells that straddle proton number 114, suggests that the Z — 114 
shell gap is not large. 

J. Concluding remarks to Section VI 

In order to judge the rehabihty of the energies of single-particle states predicted for su- 
perheavy nuclei by self-consistent mean-field theories, it is necessary to check the theoretical 
energies against the experimental ones in the heaviest nuclei where data exist. The energies 
of quasiparticle states have been calculated in a fully self-consistent manner with the CRHB 
method for ^^^Bk and 249,251q£ compared with experiment. The calculated single-particle 
energies depend on the Lagrangian parameterization; NLl, NL3 and NL-Z provide good de- 
scriptions of the measured energies, whereas NLSH and NL-RAl do not. For the former 
set, the quasiparticle energies are generally reproduced for most orbitals within ^ 500 keV. 
However, for some orbitals originating from a few specific spherical subshells, the discrep- 
ancy between theory and experiment can reach 1 MeV. Empirical shifts of the energies of 
these orbitals can be introduced to fit the experimental data. Including these shifts, the next 
spherical shell gaps beyond ^°^Pb are predicted ai, Z — 120 and N — 172, 184; no gap is seen 
ai Z = 114. The occurrence of some sizeable discrepancies in single-particle energies calls for 
an improved Lagrangian parametrization, which can better describe single-particle energies 
and, thereby, give more reliable predictions about the properties of superheavy nuclei. 

VII. CONCLUSIONS 

The cranked relativistic Hartree-|-Bogoliubov theory has been apphed for a systematic 
study of the nuclei around ^^^No, the heaviest elements for which detailed spectroscopic data 
are available. The deformations, rotational response, pair correlations, quasiparticle spectra, 
shell structure and the two-nucleon separation energies have been studied. The goal was to 
see how well the theory describes the experimental data and how this description depends 
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on the RMF parametrization. Although the relativistic mean field theory has been used 
extensively for the predictions of the properties of superheavy nuclei, it has not yet been 
demonstrated how well it describes spectroscopic data in the heaviest nuclei, which are the 
gateway to superheavy nuclei. The present investigation provides a basis for better judging 
the reliability of extrapolations to superheavy nuclei. 

The calculations with the NL3, NLl, NL-RAl and NL-Z parameter sets reproduce well 
the experimental quadrupole deformations of the Cm, Cf, Fm and No nuclei, whereas the 
NLSH set underestimate them. 

In order to reproduce the moments of inertia in the A ~ 250 mass region, the strength 
of the DIS Gogny force in the particle-particle channel has to be attenuated by 12%. In 
contrast, the moments of inertia of lighter nuclei can be well described with the full strength 
DIS force. 

With the attenuated DIS force, the rotational response is well described. In ^^^'^^"^No 
nuclei, the ahgnment of the proton ^13/2 (7r[633]7/2) and neutron ji5/2 (t^[734]9/2) pairs takes 
place simultaneously at ~ 0.31 MeV. While the crossing frequency depends only weakly 
on the RMF parametrization, the gain of aligned angular momentum at the band crossing 
and the sharpness of the band crossing are more sensitive to it. The moments of inertia at 
low spin in the Cm, Cf, Fm and No isotopes and their variations with nucleon number are 
reproduced. 

The two-particle separation energies are best described by the NL3, NLSH and NL-RAl 
parametrizations, which were derived by fitting experimental information on neutron-rich 
nuclei. The calculated deformed shell gaps occur at nucleon numbers which may deviate by 
as much as 4 from those observed in experiment. 

The quasiparticle-states calculated for odd-A nuclei are the same as those identified in 
experiment. For many states, the difference between experimental and theoretical energies 
calculated with the NLl and NL3 sets is less than 0.5 MeV, but may reach 1 MeV in some 
cases. The spectrum is less compressed in the calculations as compared with experiment, 
which reflects the low effective mass of the RMF theory. Inchision of particle-vibration 
coupling may correct that. The agreement between experiment and theory can be considered 
quite good, considering that the 6 or 7 parameters of the RMF theory have been adjusted to 
the ground-state properties of spherical nuclei, without taking into account the experimental 
information on the single-particle states. 

Concerning the predictions for superheavy nuclei we conclude the following. 

(i) Among the investigated RMF sets, NLl, NL3 and NL-Z provide best description 
of single-particle states so they seem to be most promising for the study of superheavy 
nuclei. The corresponding self-consistent calculations predict as likely candidates for magic 
numbers = 172 and A^ = 184 for neutrons and Z = 120 for protons. No significant shell 
gap is found for Z = 114. These conclusions take into account the possible shifts of spherical 
subshells that are suggested by the discrepancies between calculations and experiment found 
in our analysis of deformed odd- mass actinide nuclei. 

(n) NL-SH and NL-RAl, which are the only RMF sets predicting Z — 114 as a magic 
proton number, provide poor descriptions of single-particle states and thus are not considered 
as reliable for study of superheavy nuclei. 

(iii) Experimental studies of deformed odd nuclei with proton and neutron numbers up 
to Z 108 and A^ 164 may lead to observation of the deformed states emerging from 
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the high-j spherical sTibshclls located above Z — 120 and N — 184. Their observation will 
provide a crucial constraint on the magic numbers. 

(iv) The study of deformed states will not provide access to a number of low-j subshells, 
which largely define whether Z — 120 Z — 126 and N — 172 or N — 184 are magic 
numbers. 

(v) More systematic studies of the splitting between the 7r[521]l/2 and 7r[521]3/2 de- 
formed states, which originate from the 7r2/5/2 and 7r2/7/2 spherical subshells, may provide 
more stringent information on whether a shell gap exists at Z = 114. 

The present results demonstrate the hmitations of adjusting the RMF parameters only to 
the bulk properties of spherical nuclei and may point to missing components in the effective 
Lagrangian. A new fit to both the bulk and single-particle properties should lead to a more 
accurate theory. 

Like many Hartree(-Fock) calculations based on effective interactions, the RMF theory 
underestimates the single-particle level density. This indicates that some important mecha- 
nism is missing, which may be the particle- vibration coupling. 
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IX. APPENDIX A: THE QUANTITY 62n{Z,N). 

In order to understand the quantity S2n{Z,N) related to the second derivative of the 
binding energy as a function of nucleon number better, we first discuss the case when pairing 
is neglected. Then we perform a detailed analysis with pairing included in the CRHB-I-LN 
framework. 

Figure 29 compares 52n{Z, N) obtained in the RMF calculations without pairing with 
2Esp-GAP, where Esp-gap is the energy gap between the last occupied and the first un- 
occupied single-particle level in the {Z,N) system. One can see that 62n{Z,N) is shifted 
down by 0.561qjo MeV with respect to 2Esp-gap- This shift can be understood within 
the shell-correction method [63,102,103], in which the total energy of the system in the 
absence of pairing is given by 

Etot = ELD + E:f^ + E:^ (18) 

where Eld is a liquid drop energy and Egh is a shell energy (superscript v stands for neutrons, 
TT for protons). The latter is given by 

Esh ^ 2 ^ a - 2 j\g{e)de ^ 2 J2 - E (19) 

i—occ i—occ 
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where e is a single-particle energy, g{e) the smeared level density and E the Strutinsky 
smoothed sum. In this equation A is the Fermi energy corresponding to g{e) and is deter- 
mined from the condition of number conservation 



rX 

N^2 



I J ~g{e)de (20) 
Neglecting the variations in when the neutron number changes, one can write 

52n{Z, N) = 2ESP-GAP - SliZ, N) + Sl^'^ (Z, N). (21) 

Eld{Z, N) and E{Z, N) are smooth functions, which weakly depend on the particle number 
[63], and the 5^°(Z, iV) and 62n{Z,N) quantities related to their second derivatives as 
a function of nucleon number are nearly constant within the considered interval. Hence, 
d2n{Z, N) differs from 2Esp^gap by nearly constant. Although the shell-correction method 
is an approximation to the fully variational many-body approach such as the RMF theory, 
it elucidates the main physics in a simple way. This example clearly shows that 52n(-^, N) 
is not a direct measure of 2Esp-gap- 

The pairing smoothes the variations of S2n{Z,N) because there is gradual change of 
occupation numbers from 1 to 0. Comparing Figs. 29 and 30 one sees that the pairing 
reduces the height of the maximum of 62n{Z, N) at the N — 148 shell gap by approximately 
a factor of two and increases S2n{Z, N) at the neutron numbers away from it. This illustrates 
that the values of S2n{Z, N) cannot be taken as direct measure of 2Esp-gap since they are 
strongly dependent on pairing; see also Ref. [104] and Fig. 17b for a comparison of the results 
of CRHB and CRHB+LN calculations. 

Let us now consider the chain of the Fm isotopes within the CRHB-I-LN theory as an 
illustrative example how the 62n{Zi N) is built from different contributions. If we neglect 
the spurious ccnter-of-mass correction, which in harmonic oscillator approximation docs not 
contribute to 52n{,Z, N), then the total energy in the laboratory frame is given by (see Eqs. 
(21-24,43) in Ref. [30] for details) 

-\ga Idra{r)ps{r) - \g^ j drujQ{r)p^{r) } = Eg+v 



\gddr \lg^a\r) + \g^a\r) 



— E^NL 



- ypJdrpo{r)pi^{r) \ ^ E, 

- I e / drAo{r)ff,{r) \ = Ecoui (22) 
+ Tr{hDp) ]^Epart 

-X2{(ANy) }^Eln 

where first four terms represent the contributions from bosonic degrees of freedom, while last 
three terms from fermionic degrees of freedom. The Es+v is the sum of the energies of the 
fields associated with the linear part of the cr-meson and the a;-meson. This sum represents 
the main part of the nucleonic potential [23,32]. The E^nl term is the energy of the non- 
linear part of the cr-meson, while Ep is the energy of the p-meson field and the Ecoui is the 
energy of the Coulomb field. Finally, the Eparti Epairing, E^j^ terms are the particle and the 
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pairing energies as well as the energy correction entering into particle-number projection by 
means of the Lipkin-Nogami method, respectively. 

Based on Eq. (22), one can express S2n{Z, N) as a sum of the contributions of different 
terms of the RMF Lagrangian 

5,^{Z, N) = <5f+^(Z, N) + ^(Z, N) + 6L{Z, N) + 6g"-^{Z, N) 

+ 5Z\Z, N) + N) + 5^^{Z, N). (23) 

An analysis of these contributions to 52n{Z, N) is presented for the chain of the Fm iso- 
topes in Fig. 30. The largest contributions to S2niZ,N) come from the particle energies 
i^2n'^iZ, N)) and from the main part of the nucleonic potential (5f^^(Z, A^)). They are 
generally in opposite phase as a function of neutron number and thus they cancel each other 
to a large extent. It is interesting to see that the maximum (in absolute value) of these 
contributions is located at neutron number N — 152, while the large shell gap is seen in 
the single-neutron spectra at A" = 148 (see Fig. 19). It is difficult to understand why the 
maximum of S2n ^{Z, N) and 5f,^^(Z, A^) does not correlate with the A^ = 148 shell gap, 
but a plausible reason is related to the trend of deformation changes. The calculated /?2- 
deformation increases at neutron number A^ = 138 — 148, and then decreases at A^ > 150, 
see Fig. 14a. 

The contributions to the S2n{Z, N) coming from the p-meson {52n{Z, N)), the non-linear 
self-coupling of the a-meson {62n'^{Z,N)), the Coulomb potential (5^°"'(Z, A^)) and the 
pairing interaction (52^*""^(Z, A^)) are non-negligible and at some particle numbers some of 
them are comparable with the size of the total S2n{Z, N). 

It is interesting that total S2niZ, N) can become negative, as seen in the Fm isotopes at 
A^ = 168 in the CRHB-I-LN calculations with the NL3 and NL-RAl parametrizations; see 
Figs. 30b and f. This is a region of neutron numbers where the deformation changes are 
considerable; see Fig. 14. This result reflects their importance in the definition of 62n{Z,N) 
and again underlines the fact that many factors beyond the size of the shell gap contribute 
to 62n{Z, N). Since by definition the shell gap has to have positive value and because of the 
reasons mentioned above, it is clear that S2niZ, N) cannot be a direct measure of the shell 
gap. 

However, this quantity, being related to the second derivative of the binding energy as 
a function of nucleon number, is more sensitive to the local decrease in the single-particle 
density associated with a shell gap than the two-nucleon separation energy S2n{Z, N). 
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FIG. 2. Single-proton energies in ^ No obtained at the equilibrium deformation in the 
CRHB+LN calculations with different RMF parametrizations. Solid and dashed lines are used for 
positive and negative parity orbitals, respectively. The Xln values are shown by the long-dashed 
line with solid circles. The single-particle orbitals are labeled by means of the asymptotic quantum 
numbers [NnzA.]U (Nilsson quantum numbers) of the dominant component of the wave function. 
The asterisk (*) at the Nilsson label indicates that the wave function is fragmented and the weight 
of dominant component is below 50%. In this case, the Nilsson label does not characterize the wave 
function but is an indicator of the position of the orbital within the [N]Q, group (see Ref. [63]), 
where N is the main oscillator quantum number of the dominant iV-shell and Q the projection of 
total angular momentum onto the symmetry axis. N can be considered as a good approximate 
quantum number in almost all cases, since typically the weight of a specific A^-shcll in the structure 
of the wave function is at least 85% or larger. However, the 7r[402]3/2 and 7r[651]3/2* orbitals are 
strongly mixed by AA'^ = ±2 interaction. 
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FIG. 3. The same as Fig. 2 but for the single-neutron energies. 
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j(^) and values, respectively. In the calculations, the DIS Gogny force is attenuated by the 
scaling factor / given in Table I. The results of the calculations with NL3 and NLl parametriza- 
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FIG. 5. Proton (top panels) and neutron (bottom panels) quasiparticle energies corresponding 
to the lowest configuration in ^^^No. The CRHB+LN calculations have been performed with the 
NLl (left panels) and NL3 (right panels) parametrizations. The letters 'p' and 'h' before the 
Nilsson labels are used to indicate whether a given Routhian is of particle or hole type. 
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FIG. 6. Aligned angular momenta (in units h) versus rotational frequency Qx in ^^^No. Pro- 
ton, neutron and total {Jx) are shown. Small solid circles on lines show the frequencies at 
which the CRHB+LN calculations have been performed. The experimental quantities, defined 
as < Jj; >= + 1) ~ / + 1/2 (see Ref. [66]), are shown by open large circles. The crossing 

frequency Qx ~ 0.32 MeV is indicated by a thin vertical line. 
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FIG. 7. The results of CRHB+LN calculations with the NL3 parametrization for '^^^No. Mass 
quadrupole moments Qq, mass hexadecapole moments Q40, and 7-deformation as a function of 
rotational frequency fix are given on panels (a),(b) and (c), respectively. Proton and neutron 
pairing energies are shown in panel (d). 
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FIG. 8. Kinematic moments of inertia J^^-* (total, neutron and proton contributions) in the No 
isotopes as a function of the neutron number N. The CRHB+LN calculations are carried out with 
the NL3 parametrization at = 0.02 MeV. The vertical line shows the position of the = 148 
shell gap. 
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FIG. 9. Proton and neutron pairing energies Epairing = ~2^''(^^) obtained in the CRHB+LN 
calculations with the NLl and NL3 parametrizations for the Fm {Z = 100) isotopes and the 
= 152 isotones. 



42 




) I I I I I I I I I I I I I I I 1^ 

140 146 152 158 164 170 

Neutron number N 



FIG. 10. Kinematic moments of inertia J^^^ in the Fm isotopes as a function of the neutron 
number A^. The results of the CRHB+LN calculations at = 0.02 MeV with different RMF 
parametrizations are presented by lines. The results with the NL-Z parametrization follow those 
with NLl, but are systematically lower by ~ 1 MeV~^. The results of the calculations of Ref. 
[68] within the macroscopic+microscopic (mac+mic) method are shown by solid lines with open 
circles. Experimental data are shown by solid circles. 
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FIG. 11. The same as in Fig. 10, but for Cm, Cf and No isotopes. The results with NL-RAl 
almost coincide with those for NL3 that are displayed. 
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FIG. 12. The same as Fig. 10, but for the N = 152 isotones as a function of proton number Z. 
The results with NL-RAl are very close to those with NL3. 
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FIG. 13. The calculated (lines) and experimental (circles) deformation parameters [52 (top 
panel) and calculated mass hexadecapole moments Q^q (bottom panel) in the chain of iV = 152 
isotones. The experimental values of [32 obtained in the direct measurements [73] are shown by 
solid circles, while those deduced from the 2+ — > 0^ transition energies, with the prescription of 
Ref. [75], are given by open circles. Since the results with NL-RAl for Q4Q coincide with those 
with NL3, they are not shown in the bottom panel. 
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FIG. 14. The same as in Fig. 13, but for the chain of the Fm isotopes. The values of P2 
calculated with NLl and NL-Z are very close to each other, thus the values obtained with NL-Z 
are omitted. The results for (^40 obtained with NL-RAl and NL-Z are very close to those with 
NL3 and NLl. 
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FIG. 15. The same as in Fig. 13, but for the chains of the Cm, Cf and No isotopes. 
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FIG. 16. The two-neutron separation energies S2n{Z,N) obtained in the CRHB+LN calcula- 
tions for Fm {Z = 100) isotopes with different RMF parametrizations. Solid circles are used for 
experimental data, while open symbols for the theoretical results. 
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FIG. 17. The quantity 62n{Z,N) for Fm nuclei. The experimental data (solid circles) are com- 
pared with the results (open symbols) obtained in the CRHB+LN calculations with the indicated 
RMF parametrizations in panels (a,c-g). The results of the calculations without LN are shown in 
panel (b) by the dashed line. The experimental error bars are shown in panel (g). 
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FIG. 18. Experimental 52n{Z,N) values, shown by open symbols for Cm (Z=96) (panel (a)) , 
Cf (Z=98) (panel(b)), and No (Z=102) (panel(c)) nuclei. In all panels the experimental values for 
Fm (Z=100) nuclei are shown by solid circles in order to indicate the variations of 52n{Z,N) with 
change of proton number Z. 
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FIG. 19. Single-neutron energies in the Fm isotope chain as a function of neutron number N 
obtained at the equihbrium deformation in the CRHB+LN calculations with the NL3 parametriza- 
tion. Solid and dashed lines are used for positive and negative parity orbitals, respectively. The 
Xln values are shown by solid line with solid circles. For other details, see caption of Fig. 2. 
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FIG. 20. The same as in Fig. 16, but for the two-proton separation energies S2p{Z, N) obtained 
for the = 152 isotones. 



53 



N=152 isotones 




1.0 



V- V NL3 

□ u NL1 

_J I ^ \ \ I L 



A- — -A NLSH 
0-—0 NL-RA1 



1 



1 



1 



92 96 100 104 108 96 

Proton number Z 



100 104 108 



FIG. 21. The quantity 62p{Z, N) for the chain of iV = 152 isotones obtained in the CRHB+LN 
calculations with indicated RMF parametrizations. Solid circles are used for experimental data, 
while open symbols for theoretical results. The experimental error bars are shown in panel (b). 
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FIG. 22. The same as Fig. 19, but for the single-proton energies. 



55 



2000 



1500 - 



CO 1000 
g 

CD 

500 h 



- 



-500 - 



-1000 

CD 

-1500 I- 



-2000 



NL3+LN NL3 



EXP. NL1 



[752]3/2* 
[725]11/2 



[62011/2 \ [61519/2 

[761]1/2* 



[622]3/2 



[615]9/2 




[622]5/2 



[624]7/2 



[631]1/2 



249 



Cf 



[725]11/2 



[615]9/2 
[761]1/2* J 



[613]7/2 [620]1/2 
[734]9/2 [622]3/2 




[743]7/2 



[501]1/2 



FIG. 23. Experimental and theoretical quasiparticle energies of neutron states in ^^^Cf. Positive 
and negative energies are used for particle and hole states, respectively. The experimental data 
are taken from Ref. [81]. Solid and dashed lines are used for positive and negative parity states, 
respectively. The symbols 'NL3' and 'NLl' indicate the RMF parametrization. The CRHB results 
shown below them were obtained with original DIS Gogny force (/ = 1.0) used in pairing channel 
and without particle number projection. 'NL3+LN' indicates results with the LN method, the NL3 
parametrization and scaling / of the strength of DIS force (given in Table I). In each calculational 
scheme, attempts were made to obtain solutions for every state shown in figure. The absence of a 
state indicates that convergence was not reached. 
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FIG. 24. The same as in Fig. 23 but for neutron states in ^^^Cf. The experimental data are 
from Refs. [82,83,89]. 
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FIG. 25. The same as in Fig. 23 but for proton states in ^^^Bk. The experimental data are 
from Ref. [80]. 
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FIG. 26. Neutron and proton single-particle energies in No. The columns marked by 'NLl' 
show the original spectra obtained with the NLl parametrization (see Figs. 2 and 3). The columns 
'NLlcor' show how the spectra are modified if the energies were shifted as discussed in Sect. VIE 
and VI D. Solid and dashed lines are used for positive and negative parity states. Deformed gaps 
are indicated. 
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FIG. 27. Proton single-particle states in a 172I2O nucleus. Columns 'NL3' and 'NLl' show the 
states obtained in the RMF calculations at spherical shape with the indicated parametrizations. 
The energy of the 1^13/2 state in the NLl parametrization is set to be equal to that in NL3, which 
means that the energies of all states in NLl (last column) are increased by 0.78 MeV. The columns 
'NLScor' and 'NLlcor' show how the spectra are modified if empirical shifts were introduced based 
on discrepancies between calculations and experiment for quasiparticle spectra in deformed ^^^Bk 
(see Sect. VIE for numerical values). Solid and dashed lines are used for positive and negative 
parity states. Spherical gaps at Z = 114 and Z = 120 are indicated. 
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FIG. 28. The same as Fig. 27, but for neutron single-particle states. The energies of all states 
obtained with the NLl parametrization (last column) are increased by 0.76 MeV in order to have 
the same energies of the 2(79/2 states in the second and third columns. Spherical gaps at = 172 
and N = 184 are indicated. 
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FIG. 29. Comparison of the 52n{Z,N) quantity with the twice the size of the shell gap 
2Esp-GAP, which is the distance between the last occupied and first unoccupied orbitals. Both 
quantities are calculated with the NL3 parametrization and without pairing. 
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FIG. 30. The contributions 4n ^) of the different terms of the CRHB+LN theory to the 
quantity 52n{Z, N) as a function of the neutron number for the chain of the Fm isotopes. The 
contribution is not shown since it typicahy Ues in the range from —30 keV up to +30 keV. 
The only exceptions are = 148 and = 162, where = 83 and 62 keV, respectively. 
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TABLES 

TABLE I. The scaling factors / of the Gogny DIS force (see Eq. (10)) used for different 
parametrizations of the RMF Lagrangian in the CRHB+LN calculations. 



Parametrization 


NLl 


NL-Z 


NL3 


NL-RAl 


NLSH 


/ 


0.893 


0.880 


0.864 


0.861 


0.876 



TABLE IL Three-point indicators of the odd-even staggering of binding energies 
A(3)(iV) = [B{N - 1, Z) + B{N + 1, Z) - 2B{N, Z)], where B{N, Z) is the (negative) bind- 

ing energy of a system with N neutrons and Z protons. An analogous proton indicator /S.^'^\Z) 
is obtained by fixing the neutron number N and replacing N by Z. Column 1 indicates the type 
of indicator (proton /s!^'^\Z) or neutron /S.^^\N)) and nucleus with Z (proton indicator) and 
(neutron indicator). Columns 3 and 4 give the results obtained in CRHB calculations with / = 1.0 
for the DIS force and with the NL3 and NLl parametrizations, while columns 5 and 6 give those 
obtained in CRHB-I-LN calculations with / values given in Table L 



A(3)(...) 


exp 


NL3 


NLl 


NL3-KLN 


NLl-FLN 


1 


2 


3 


4 


5 


6 


A(3)(Z) p^'^Bk] 


0.399 


0.516 


0.515 






A(3)(Ar) [249cf] 


0.519 


0.481 


0.559 


0.458 


0.515 


A(3)(Ar) [25icf] 


0.531 


0.491 


0.605 
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TABLE III. Spherical subshells active in superheavy nuclei and their deformed counterparts 
active in the A ~ 250 mass region. The left column shows the spherical subshells active in the 
vicinity of the "magic" spherical gaps (Z = 120, N = 172). Their ordering is given according to 
the RMF calculations with the NL3 parametrizations in the 172 120 system (see Figs. 27 and 28). 
Although the gaps depend on the specific RMF parametrization, the same set of spherical subshells 
is active with other parametrizations (see, for example, Fig. 4 in Rcf. [17]). The right column shows 
the deformed quasiparticle states observed in ^97Bki52 [80] and ^gg'^^^Cfisi^isa [81-83]. The bold 
style is used for the states which may be observed in nuclides with N ~ 162 and/or Z ^ 108. The 
symbols 'N/A' (not accessible) are for the deformed states which typically increase their energy 
with deformation and thus are not likely to be seen experimentally. 



Spherical subshell 


Deformed state 


Proton states 




-1 ' 

7rljl5/2 


TT [77UJ1/2 


7r3pi/2 


IN /A 


7r3p3/2 


AT / A 
IN /A 


7rl?;ii/2 


TT [651] 1/2 


rz 1 on 

Z = izU 




7r2/5/2 


fCTOlll /O 

7r[521Jl/2 


7r2/7/2 


7r[521J3/2, 7r[530Jl/2 


7rlZiq/9 


TT [642] 5/2, TT [633] 7/2, tt [624] 9/2 


7r3si/2 


TT [400] 1/2 


nlhg/2 


7r[514]7/2 


Neutron states 




iylkn/2 


V [880] 1/2 


iy2hu/2 


V [750] 1/2 


^^1^13/2 


i/[761]l/2 


N = 184 




iy4si/2 


N/A 


vM^l2 


[620] 1/2 


vM^l2 


N/A 


N = 172 




I^2fif7/2 


zv[622]3/2 


^"^99/ 2 


z/[622]5/2, zy[613]7/2, zy[604]9/2 


l^ljl5/2 


z/[734]9/2, z/[725]ll/2 


Z/lZll/2 


z/[615]9/2, z/[624]7/2 
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TABLE IV. Additional binding (AE^tvm) induced by nuclear magnetism for different neutron 
(^^^Cf) and proton (^"^^Bk) quasiparticlc states obtained in the CRHB calculations with the NL3 
parametrization and full strength (scaling factor / = 1.0) of the Gogny force. The quantity AEnm 
is defined as the difference of binding energies obtained in the calculations with and without nuclear 
magnetism. As tested for a number of states, the change of the scaling factor / to the one given 
in Table I and/or the use of the LN method modifies AEnm only marginally. 



neutron state 


AEnm (keV) 


proton state 


AEnm (keV) 


[734] 9/2 


-36 


7r[521]l/2 


-16 


z/[615]9/2 


-55 


7r[514]7/2 


-35 


i/[624]7/2 


-56 


7r[633]7/2 


-22 


z^[622]3/2 


-27 


7r[624]9/2 


-23 


z^[622]5/2 


-33 


7r[521]3/2 


-27 


[734] 7/2 


-37 


7r[523]5/2 


-33 


z^[613]7/2 


-29 


TT [642] 5/2 


-23 


[725] 11/2 


-34 






z/[761]l/2* 


-69 






i/[752]3/2* 


-53 
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